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ABSTRACT 

We present Very Large Array 21-cm observations of the massive edge-on spiral 
galaxy NGC 5746. This galaxy has recently been reported to have a luminous 
X-ray halo, which has been taken as evidence of residual hot gas as predicted 
in galaxy formation scenarios. Such models also predict that some of this gas 
should undergo thermal instabilities, leading to a population of warm clouds 
falling onto the disk. If so, then one might expect to find a vertically extended 
neutral layer. We detect a substantial high-latitude component, but conclude 
that almost all of its mass of 1.2 — 1.6 x 10^ Mq most likely resides in a warp. 
Four features far from the plane containing about 10^ are found at velocities 
distinct from this warp. These clouds may be associated with the expected infall, 
although an origin in a disk-halo flow cannot be ruled out, except for one feature 
which is counter-rotating. The warp itself may be a result of infall according to 
recent models. But clearly this galaxy lacks a massive, lagging neutral halo as 
found in NGC 891. The disk HI is concentrated into two rings of radii 1.5 and 
3 arcminutes. Radial inflow is found in the disk, probably due to the bar in this 
galaxy. A nearby member of this galaxy group, NGC 5740, is also detected. It 
shows a prominent one-sided extension which may be the result of ram pressure 
stripping. 
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Subject headings: galaxies: ISM — galaxies: spiral — galaxies: kinematics and 
dynamics — galaxies: formation — galaxies: evolution — galaxies: individual 
(NGC 5746, NGC 5740) 



Introduction 



Gaseous thick disks or halos in spiral galaxies hold promise for answering key ques- 
tions about galaxy formation and evolution. These vertically extended layers are multi- 



phase, wi th detailed s tudies of the nearest edge-on galaxies revealing halo s of neutral 
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Especially for the extended DIG, X-ray, dust and radio continuum components, it has 
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origin of these layers in a star-formation-dri ven disk-halo cycle, which has been th eoretically 



modeled as a general galactic fountain flow (IShapiro fc Fieldlll976l : lBregmanlll980l ) with later 



models incorporating the fact that mass and energy input into t he halo may efficiently oc cur 



through localized structures such as supershells and chimneys ( iNorman fc Ikeuchilll989l ). 



As for the Milky Way, the WHAM survey (IHaffner et al.l l2003r) has characte rized in 
detail the so-called Reynolds layer of vertically- extended DIG (IRevnolds et al.lll973r) . reveal- 
ing further instances of possible superbubbles ([Reynolds et al.ll200ll ). Very relevant to this 
paper are the High Velocity Clouds (HVCs) and Intermediate Velocity Clouds (IVCs). For 
the latter, a few are known to be 0.3-4 kpc above the midplane, with met ali i cities close to 



solar, and it is quite possible that these clouds originate in a disk-halo flow (IWakkerl |2004J . 
and references therein). The situation for the HVCs is different, h owever. Not considering 
the contribution from the Magellanic Stream (jPutman et al.ll2003l ). some well studied com- 
plexes are found to be many kpc from the plane (e.g. 8-10 kpc for part of complex A, >4 
kpc for complex C; IWakkeii l2004j and references therein). Although information is scarce. 
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their metallicities are also well below solar (e.g. ICoUins et al.l 120071 ). suggesting that such 
clouds are not part of a disk-halo flow but may be extragalactic clouds that are mixing with 
metal-enriched gas in the halo as they fall onto the disk. 

Such a population of infalling warm clouds is a prediction of recent models of galaxy 
formation, in which galaxies are still growing by this mechanism. The idea that g alaxie s form 



by gas cooling out of a shock-heated hot halo was presented by IWhite fc Rees 



hot gas should be therma ll y unstable and giv e n to f ragmentation (IField 



1965 



1978). This 



Fall fc Rees 



19851 : [Murray fc LinI Il990l ) . iMaller fc BuUockl (120041 ) explored the consequences of such an 
inflow of fragmenting material. The remaining hot gas has such a low density that it can 
support itself against infall for a long time. This may explain the "over-cooling" problem 
in galaxy formation models that leads to overly massive galaxies and the find ing that most 



of the baryons in the universe never ended up in galaxies (IBenson et al 



20031 and references 



ther ein). Instabilities in c ooling infiows have also been studied by iKaufmann et al.l (120061 ) 
and ISommer-LarsenI (120061 ). The cloud populations and distributions in these high- resolution 
models vary, with clouds having an uncertain but probably high ionization fraction. 

An alternative explanation to the "over-cooling" problem in the most massive galaxies 
is feedback. Reheating of the halo gas by core-collapse supernovae o r a (possibly recurrent) 
AGN prevents the gas fror a cooling ari d fiowing in (iBinneyl 120041 ). Heating by Type la 
supernovae is considered by IWang (120051 ). It is also argued that, except in the most naassive 
halos, infalling gas never heats to the virial temperature but flows in cold (iBinnevI 119771 : 
Birnboim fc Dekellbood : lBinnevl[20oi : iKeres et al.lbood : lOekel fc Birnboimll2006h . 



The rotation of gaseous halos may also provide clues as to their origin and the physical 
processes occurring within them. In recent years, the manner in w hich rotation, speed s c hange 
as a function of heigh t has begun to be characterized in the DIG dHeald et al.ll2006al lbl. I2OO7I : 
Oosterloo et al.l 120071 ). Heald and coworkers have measured the gradient in rotation speed 
with height [dVrot/dz) in three edge-ons which form a decreasing sequence of star forming 
activity and DIG scale-height: NGC 5775 (-8 km s^^ kpc^^), NGC 891 (-17 km s'^ kpc^^), 
and NGC 4302 (-30 km s'^ kpc'^ 

DIG scale-height, show much less range: -15 to - 25 km s~^ (scale he ight) 



The gradients, when expressed in terms of km s ^ per 



The authors 

show that simple ballistic models of disk-halo flow (jCoUins et al.l 120021 ) greatly underpredict 
the magnitude of these gradients and also predict the wrong trend with scale-height. Two 
possible resolutions are: 1) the physics of disk-halo flows are not well describe d by ballistic 
models and that hydro dynamical effect s such as pressure (IBarnabe et al.l 120061 ) or magnetic 
or viscous forces (jBenjaminl |2002| ) dominate the rotation, or 2) part or all of the gas does 
not originate in a flow from the high-angular-momentum disk, but rather from infalling low- 
angular-momentum primordial gas. For the DIG, the latter explanation must also account 
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for morphological and other connections with ongoing star formation discussed above. 



For NGC 891, dV mt/dz has also been measured for the HI and is -15 km s~^ kpc~^ 
(lOosterloo et al.ll2007l ). in good agreement with the DIG value . The authors conclude from 
this gradient that at least sort ie of the 1.2 x 10^ M^r^ HI halo must be accreted, low angular 
momentum material (see also iFraternali fc Binneyl |2006| ) . 



Furthering our understanding of both galactic infall and disk-halo cycling, as well as 
possible interactions between the two, would be made easier if one could isolate galaxies 
where one or the other origin is expected to obtain. High star-forming galaxies like NGC 
5775 presumably have gaseous halos dominated by disk-halo flows. The recent discovery of a 
bright ( Lrc = 7.3 ±3.9 x 10^^ erg s~^) X-ray halo around the massive, nearby ( 29.4 Mpc is the 



commonly adopted distance) , low star forming edge-on Sb galaxy NGC 5746 (jPedersen et al. 



20061 : iRasmussen et al.ll2006l ) pr esents a challenge to disk-halo models as the X-r ay luminosity 
clearly exceeds that expect ed (ITuUmann et al.l l2006al : [Strickland et al.l 120041 ) for a galaxy 



with no detected DIG halo (jRandlll996l: in an imag e with an Emission Measure noise level of 



3.7 pc cm~^) and little star formation. iRandl (119961 ) roughly characterized the star formation 
rate per unit disk area of many edge-ons using the far infrared luminosity measured by the 
Infrared Astronomical Satellite (IRAS) divided by the optical disk area: Lfir/D\^ (values 
calculated from the NASA/IPAC Extragalactic Database!) . NGC 5746 is one of the lower of 
the 16 edge-ons thus characterized, at 4 x 10^^ erg s~^ kpc~^. Rather, the X-ray luminosity 
puts NG C 5746 on the ex pected steep relationship oc VJg^ for hot halos in galaxy formation 



models (IToft et al.ll2002l ). and, in such an interpretation, it is only because of the very high 
rotation speed (measured in this paper to be about 310 km s~^) that this residual hot halo 
can be detected at all. The metallicity of the hot gas is found to be low, at about 0.0 4 solar 
(|Rasmussen et al.ll2006l ). but is uncertain and could be biased toward low values (see iBuote 
2000l ). But overall, this galaxy may well be an attractive test case where disk-halo cycling 
is minimized. 

There are some cavea ts to this interpretat ion, however. First, NGC 5746 is in a group 
of 26 cataloged members (iGiuricin et al.ll2000l ). and there may well have been interactions 
in the past leading to gas in the group environment. Second, although there is little ongoing 
star formation, it may be that Type la supernovae or an AGN create a ho t wind or provide 



energy to keep a pre-existing halo heated, as discussed above. However, IRasmussen et al. 



(120061 ) argue that the likely Type la SN rate in a galaxy like NGC 5746 would be far too low to 
explain the X-ray luminosity. They also argue that there is little evidence for at least current 



^The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, Cali- 
fornia Institute of Technology, under contract with the National Aeronautics and Space Administration. 
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nuclear activity in this galaxy, and that a soft, thermal X-ray halo would not be expected 
for an AGN outflow aiiy way. We do note, however, that NGC 5746 is classified as a LINER 
by ICarrillo et al.l (119991 ). and there is a compact X-ray source (iGonzalez-Martm et al.ll2006l ) 
at the center, indicating at least low level a ctivity at present, although no compact radio 
source has been detected ( iNagar et al.ll2005l ). If the metallicity of the hot gas is indeed as 
low as 0.04, these alternative sources are not very likely. 

If this X-ray halo is a residual of galaxy formation and thus an indication of missing 
baryons, then an obvious question is whether the predicted thermal instabilities are occur- 
ring, leading to an infalling warm component. If it exists, the lack of a DIG halo may simply 
mean that there are few disk sources capable of significant ionization. But there may be a 
neutral halo. Either a detection or an upper limit will constrain galaxy formation models 
and the "missing baryon" question, and shed light on the origin of the Milky Way's HVCs, 
more so if individual clouds can be seen. A detection of a significant neutral halo would also 
be challenging to explain in terms of a disk-halo cycle origin. 

We have therefore observed NGC 5746 in 21-cm emission with the Very Large Array 
(VLA), as described in §2. We analyze high-latitude emission in §3, and end with a brief 
discussion in §4 in terms of the theories of the halo gas discussed above. 



2. Observations 

NGC 5746 was observed in the C array of the VLA on 2007 January 7, 8, 12 and 13. 
Phase calibration was achieved through observations of VLA calibrator 1445-1-099 about 
every 30 minutes. Observations of 3C286 and 3C48 were used for flux and bandpass calibra- 
tion. Sixty three spectral channels were employed, centered at 1725 km s~^, with channel 
width 20.85 km s~^, while online Hanning smoothing yielding a velocity resolution of twice 
this width. A total of about 28.3 hours were spent observing NGC 5746. Eight of the 27 
antennae were unavailable for most of the final track, while two were absent on the other 
three dates. Smaller amounts of data were also lost due to high winds, interference, and 
equipment failures. Data were inspected for high amplitudes and any suspect data excised. 
Some small baseline corrections for a few antennae were made, based on later observations 
in the same configuration. The calibration of all four tracks is of very high quality. 

High amplitude visibilities were clipped before mapping. Dirty maps of each track were 
made to check for additional problems before concatenating the four uv datasets. Continuum 
was subtracted in the wf-plane using line-free channels at either end of the spectrometer with 
the AIPS tasks UVSUB and UVLSF. This was very successful. The AIPS task IMAGR was 
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used to produce clean maps using clean boxes that covered the emission from NGC 5746 as 
well as NGC 5740, which is also detected, albeit somewhat outside the 31.9' FWHM primary 
beam. The pointing center is R.A. 14^ 44°^ 56.4% Decl. 1° 57' 16" (J2000). The primary 
data cube discussed here was made with no uf-tapering and has a resolution of 15.5"xl4.8" 
(2.2 X 2.1 kpc) at a P.A. of -3.6°. A cube with 61.3"x58.1" resolution (8.7 x 8.3 kpc), at 
a P.A. of 45.5°, created via wf-tapering is also discussed here. All maps have 1024x1024 
pixels of 3" size, and were made with uniform wf-weighting with the IMAGR 'ROBUST' 
parameter equal to 0. The noise in a single channel of the full-resolution cube is 0.23 mJy 
(beam)"^. No primary beam correction was made for maps of NGC 5746. Primary beam 
attenuation reaches 4% at the ends of the major axis. Maps of NGC 5740 are corrected for 
primary beam attenuation; the response at the center of this galaxy has dropped to 0.36, and 
thus structure may not be accurately mapped. No other galaxies in the NGC 5746 group 
are within the mapped region. Observed velocities are on the heliocentric scale. The noise 
corresponds to an HI column density of 4.6 x 10^^ cm~^ averaged over one resolution element 
for optically thin gas. In a zeroth-moment map made using all data and all channels, a 5a 
detection of a poin t source corresponds to 2.3 x 10^ Mq, using the conversion from (e.g.) 



Zwaan et al.l (119971 ). All masses are of total atomic gas and include a correction of 1.36 for 



helium content. 

Finally, the GIPSY task BLOT was used to blank emission-free regions in each channel, 
reducing the noise in moment maps. Moment maps were made with various strategies to 
eliminate as much noise as possible while retaining as much emission as possible that appears 
to be real in the channel maps and position- velocity (pv) diagrams (given its signal-to-noise 
ratio and continuity over multiple velocity channels). The zeroth-moment map presented 
here for the full-resolution cube includes all emission in the blotted cube above the la level 
in each pixel in two consecutive channels. For the 60" resolution cube, such a strategy 
introduced too many features in the zeroth-moment map that did not appear real in the 
data cube, and thus a 2.5a cutoff was used. A concern with the velocity resolution of 42 km 
s~^ is that many real spectrally narrow features may be rejected by this strategy. 



3. Results 

Channel maps of the full- resolution cube are shown in Figure 1. Zeroth- and first- 
moment maps made from the blotted full-resolution cube are shown in Figure 2. The zeroth- 
moment map is overlaid on an Digitized Sky Survey red image in Figure 3. Also labeled 
in Figure 3 are the four extraplanar features discussed in §3.5. The emission from NGC 
5746 is dominated by what appears to be a highly inclined disk at a position angle (PA) of 
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350° (measured CCW from N to the receding side of the major axis) with a large central 
depression of about 2' diameter. In this se nse it is somewhat r eminiscent of the Infrared Space 
Observatory ISOCAM 12/im image from iBendo et al.l (120021 ) . which has the appearance of 
an inclined ring of approximate diameter 2.3' (19.7 kpc). It is possible the HI hole is larger 
but limited resolution has made it appear smaller. 
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Fig. 1. — Channel maps from the full- resolution cube for NGC 5746. Contour levels are 3, 
6, 12, 24, 36 and 48 times the 1 a noise in column density units of 4.6 x 10^^ cm~^. Each 
map is labeled with its heliocentric velocity. 
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R.A. (2000.0) R.A. (2000.0) 

Fig. 2. — Left: Zeroth-moment map of NGC 5746 from the full-resolution cube. Contour 
levels are 1, 2, 4, 8, 17 and 32 times 10^° cm~^. The beam is shown in the lower right corner. 
Right: First-moment map from the full-resolution cube. Contours run from 1500 to 2020 km 
s~^ in 40 km s~^ increments from south to north. 
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45"" 14^44"'55' 44°'50' 
R.A. (2000.0) 



Fig. 3. — Zeroth- moment map of NGC 5746 from the full- resolution cube overlaid on a 
Digitized Sky Survey red image. Contour levels are as in Figure 2 (left). The approximate 
locations of the centers of the four features discussed in 53.5 are also indicated. 
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At fainter levels, a component extended along the minor axis is clearly detected. It 
is resolved into individual clouds and filaments to some de gree. This is in contrast to the 
HI map of NGC 891 (see Figure 1 of lOosterloo et al.l 120071 ) where the halo is smoother in 
appearance despite better linear resolution (1.4 kpc vs. 2.1 kpc). Emission is detected up to 
about 70" (10 kpc) from the midplane. The northern end of the disk also suggests a warp. 

The total fiux in the map in Figure 2 is 33.7 Jy km s~^, yi elding a total at o mic m ass of 
9.4 X 10^ M0. For comparison, the total observed fiux found by lSpringob et al.l (120051 ) using 
the Arecibo telescope is 38.8 Jy km s~^ 

Figure 4 shows the zeroth- and first-moment maps from the 60" -resolution cube. The 
total mass in this map is 8.8 x 10^ Mq, slightly less than in the full-resolution map. We 
attribute this to the somewhat higher cutoff that was necessary to eliminate noisy features 
in the moment map at this resolution. However, it also indicates that there is very little 
diffuse, low surface brightness emission that is missed in the full resolution cube. The ends 
of the disk in this map also suggest a warp, but with the opposite sense of bending than is 
evident at lower radii in the full resolution map. 
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R.A. (2000.0) R.A. (2000.0) 

Fig. 4. — Left: Zeroth-moment map of NGC 5746 from the 1 '-resolution cube. Contour 
levels are 1, 2, 4, 8, 16, 32, 64 and 92 times 1.3 x 10^^ cm~^. The beam is shown in the lower 
right corner. Right: First-moment map from the 1 '-resolution cube. Contours run from 1460 
to 2020 km s~^ in 40 km s~^ increments from south to north. 
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3.1. Modeling the Atomic Mass Distribution and Rotation 

Our main goal here is to understand the origin of the faint, vertically extended com- 
ponent seen in Figure 2 through full modeling of the data cubes. Although we are able to 
quantitatively constrain many of the model parameters, we do not carry out a full search 
of parameter space, not least because of asymmetries in the data, but rather explore and 
constrain a few well motivated types of models. 

We begin by constraining the radial density profile at the midplane. As a first estimate, 
we employ the GIPSY task RADIAL, which fits radial profiles of the column density inte- 
grated vertically through the disk to averaged major axis emission profiles for highly inclined 
galaxies. The observed major axis emission profile is shown in Figure 5 as the solid line. 
Since in the below modeling we will use axisymmetric density distributions, we also show 
a profile where the north and south sides have been averaged to produce an axisymmetric 
profile (dashed line). The fitted axisymmetric profile is shown as the dotted line. The radial 
column density profile derived by RADIAL is shown in Figure 6 as the filled circles. The disk 
is dominated by two rings of radius 1.5' and 3' and width about 1', while the profile confirms 
the aforementioned central depression of radius 1'. RADIAL reproduces the observed pro- 
file in Figure 5 very well, except for a moderate overestimate of the emission in the central 
depression. This overestimate cannot be rectified by creating a more prominent central hole 
but is instead a result of projected emission from gas at larger radii; it presumably indicates 
an asymmetry in the galaxy. The other symbols in Figure 6 show the deviations from the 
basic profile featured in the models described below. 
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Major axis offset (arcmin) 



Fig. 5. — The solid line shows the observed major axis emission profile averaged over the 
minor-axis extent of the emission in NGC 5746. The dashed line shows the average of the 
north and south sides of the observed profile. The dotted line shows the fitted profile from 
the GIPSY task RADIAL. 
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Fig. 6. — Filled circles show the axisymmetric radial profile of the column densities integrated 
vertically through the disk derived by RADIAL for NGC 5746. This profile is used for the 
unwarped component of Models A and B. Open circles show the modification necessary to 
model the extended emission along the minor axis as a warp in Models A and B. Triangles 
show the slightly lower profile for the disk component in the lagging halo Models C, D, and 
E. The 'x' and the open squares show the column density of the ringlike lagging halo Models 
D and E, respectively. See text for further details of the models. 
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We use the major axis pv diagram (Figure 7) and the envelope tr a cing m ethod to 
estimate the rotation curve. The method is described in ISofue &: RubinI (1200 ll ). We use 



7] = 0.3 (see ISofue &: RubinI l200ll for the definition of this parameter). We assume the 
inchnation is close enough to 90° so that no correction for projection is necessary. As 
there is little modeled column density in the central 1' radius, the rotation curve is very 
poorly constrained there. In fact, long-slit spectra and a peanut bulge indicate that NGC 
5746 is clearly barred, with the bar elongated more across the line of sight than along it 
(IKuijken fc Merrifieldll 19951 : iBureau fc Freemarul 19991 ) . The bar manifests itself kinematically 
in a steeply rising component within major axis offsets of ±0.25' in pv diagrams of emission 
lines (indicative of X2 orbits perpendicular to the bar), which is clearly not seen in HI. For 
our rotation curve we have simply extended the value at = 1.15' inwards. Our adopted 
rotation curve is shown in Figure 8. The dynamical mass within R = 5' (42.8 kpc), given by 
M{R) = RV\R)/G, is 7.9 x 10^^ M©. 
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Fig. 7. — At the top, middle and bottom are shown major-axis position-velocity diagrams 
for NGC 5746 from the full-resolution cube, Model B and Model D respectively (see text for 
model details). Contour levels are as in Figure 1. 
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Fig. 8. — The filled circles show the rotation curve for NGC 5746 derived from envelope 
fitting. The open circles show the assumed extension at low radii. 
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We next need to constrain the inclination, i. For nearly edge-on galaxies, tilted-ring 
programs such as ROTCUR in GIPSY cannot do this accurately, especially if a warp along 
the line of sight or a lagging halo is suspected and contributes significantly to the velocity 
field. However, i can be estimated by assuming axisymmctry in the ringlike appearance of the 
emission at the center of the zeroth-moment map. Hence, wc use our adopted rotation curve 
and radial column density profile as inputs to the GIPSY program GALMOD (a program 
which allows construction of model galaxies with specified radial column density profiles, 
rotation curves, and various forms for the vertical density distribution) in order to generate 
a zeroth-moment map and match it to the data. The radial bins in all our models have a 
width of 15" (2.1 kpc). The models produced by GALMOD are convolved to the resolution 
of the data. From the measured separation of the two sides of the ring along the minor axis 
near the galaxy center in Figure 2, wc find that i must be about 86°. This value is confirmed 
by the modeling described below. ROTCUR was run to constrain other parameters, and 
indicates a PA of 350 ± 0.3°, a systemic velocity of 1733 ± 7 km s^^, and a kinematic center 
about 6'' or 850 pc west of the pointing center, at R.A. 14^^ 44°' 56.0", Decl. 1° 57' 16" 
(J2000). The gas is initially assumed to be in a single exponential layer with scale height 3" 
or about 400 pc. The true scale height of the bulk of the emission is not well constrained 
because of the resolution and could be substantially lower than 3" . 

We next show that the vertical structure cannot be fit by a single component at a given 
inclination, or by two exponential components. Figure 9 shows a minor-axis emission profile, 
averaged over 3' along the major axis (solid line). There are clearly tails at both ends of 
the profile. The dashed line shows our modeled component for an exponential layer with 3" 
scale height. It fits the bulk of the emission well but not the tails. The dotted line shows an 
exponential model with a scale height of 18". It fits the tails well but is broader than the 
bulk of the emission and washes out the splitting at low latitudes (which is due to the ringlike 
distribution). No combination of two exponentials were found to produce a reasonable fit. 
This result does not change if a gaussian or sech^ form is used. 
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Fig. 9. — The solid line in both panels shows the minor axis emission profile averaged over 
the inner 3' (25.6 kpc) of the major axis, in column density units. The dashed and dotted 
lines in the top panel show profiles for models as described in the text with exponential scale 
heights of 3" and 18", respectively. The dashed line in the bottom panel shows the profile 
for Model A, as described in the text. The profiles for the other models described are almost 
identical to that of Model A. 
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Instead we consider two possibilities to explain the tails. The first is a warp along the 
line of sight, and the second is a halo with a box-like (employed for ease of modeling although 
it is unphysical) profile in z. The model parameters are summarized in Table 1. 

The morphology of the tails turns out to be very much coupled with their kinematics, 
so we now consider these together by introducing pv diagrams parallel to the minor axis 
at various positions along the major axis which provide good leverage on the parameters of 
interest. These are shown in Figure 10, along with several models generated by GALMOD 
and variants on that program. The characteristic appearance of the disk component in these 
diagrams is one of narrow angular extent at velocities furthest from Vsys which broadens 
and (where the ringlike structure is evident) splits as velocities move towards Vsys- The 
faint high-latitude component has a characteristic appearance which will greatly constrain 
its morphology and kinematics. Although it varies significantly among cuts, it generally 
manifests itself as an extended, spectrally very narrow component at velocities closer to Vgys 
than the disk component, with a definite gradient in T^ej with h that is more evident in panels 
further from the galaxy center. Some east- west asymmetries are evident in this component. 
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Fig. 10. — Position-velocity diagrams parallel to the minor axis at various offsets along 
the major axis from the full-resolution data cube and various models described in the text. 
Contour levels are as in Figure 1. The offsets along the major axis are shown in the panels 
for the data. 
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We first consider a warp along the line of sight. In GALMOD, this is achieved by slowly 
decreasing the inclination of outer rings. However, in order to reproduce the minor axis 
emission profile, the column density in the warped rings had to be increased significantly 
from the major-axis-based starting point, as shown in Figure 6. The warp begins at a radius 
of 4' (34.2 kpc), where i = 84°, dropping to 79° by 5' (42.8 kpc) radius (Figure 11). The 
maximum displacement of the midplane of the outermost warped ring from the unwarped 
disk is about 5 kpc. A slightly larger warping is required for the SE quadrant (top three 
panels of Figure 10 at negative minor axis offsets) but is not modeled here. The resulting 
model, Model A, appears to match the observed minor-axis pv diagrams very well. 
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Fig. 11. — Filled circles show the run of inclination (left scale) for Models A and B. Open 
triangles show the run of PA (right scale) for Model B. Open circles show the inclination for 
Models C, D and E. 



- 25 - 



Model B is an attempt to incorporate the bending across the hne of sight evident in the 
zeroth-moment maps. Retaining the run of inchnation in Model A, the PA in this model 
decreases from 350° to 346° over the radius range 3' to 4', and then increases back to 5° by 
5' radius (Figure 11). The asymmetry introduced into the minor-axis pv diagrams (Figure 
10) provides a slightly better match to the data at high latitudes, although still probably 
underestimating the observed asymmetry. The zeroth-moment map of this model (Figure 
12) also exhibits a slight brightness asymmetry with respect to the minor axis at major 
axis distances less than 4' from the center and at heights 30-60" from the plane that is 
also apparent on the east side in Figure 2. The model appears to somewhat overestimate 
the bending to larger PAs evident in Figure 2 (and, of course, there is too much emission 
at the ends of the major axis because of the excess column density in the modeled warp). 
However, a model version smoothed to 60" (Figure 12, center), is a reasonable match to th e 
bending in Figure 4. We note that warps are often asymmetric (iGarcia-Ruiz et al.ll2002l ). 
so one should not expect a warp fitted to kinematic information parallel to the minor axis 
to match perfectly the morphology along the major axis. A second caveat to this exercise 
is that some of the bending along the major axis may be due to outer spiral structure seen 
not quite edge-on. In fact, the initial bending to lower PAs in the model is only present to 
match the morphology in the moment maps. The subsequent increase in PA is responsible 
for the asymmetries in the pv diagrams. As mentioned at the beginning of this section, for 
these kinds of reasons, we have not carried out a full exploration of parameter space but 
simply present these models as reasonable fits to the data. 
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Fig. 12. — Left: Zeroth-moment map of Model B, convolved to 15" resolution. Contour 
levels are as in Figure 2 (left). Center: Zeroth-moment map of Model B, convolved to 60" 
resolution. Contour levels are as in Figure 3. Right: Zeroth-moment map of Model D, 
convolved to 15" resolution. Contour levels are as in Figure 2 (left). 
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We next consider the box-like halo, extending to z — 50" (7.1 kpc). Such a component, 
when convolved with the beam, loses its sharp edges and, added in the right proportion to 
the disk component, reasonably reproduces the minor-axis emission profile in Figure 9, being 
almost indistinguishable from the warp model profile. The disk component column density 
is normalized downwards by a factor of 1.24 in these models. 

For the halo models, the gradient in T4e/ with b at high latitudes in Figure 10 implies a 
rotational velocity that decreases with z, i.e. dVrot/dz < 0. The spectral narrowness implies 
that such a halo cannot occupy a large range in radius, otherwise projection of the rotational 
velocity would cause broadening of the Vhei profiles in disagreement with the data. 

To illustrate these constraints, we first consider a model where the halo has the same 
radial density profile shape as the disk, and a vertical gradient in rotation velocity of -3.5 
km s~^ (arcsec)"^ (-25 km s^^ kpc^^). We refer to this as Model C in Figure 10. Although 
the amount of high-latitude emission is well matched to the data, there is little such emission 
seen at the first sohd contour level of this model. This is because the velocity profiles there 
are broader with lower peak intensities, extending to velocities much too far from Vsys to be 
consistent with the data, and we demonstrate this by including a low-level dashed contour. 
Clearly, this model is a poor match to the halo emission. 

The narrowness of the high-latitude velocity profiles puts a strong constraint on the 
radial range of the halo component, and leads us to consider annular distributions, beginning 
with one in which the entire halo is contained in one 15" ring. Model D in Figure 10 features 
a ring of central radius Ro = 180", with dVrot/dz — —3.5 km s~^ arcsec"^, and a necessarily 
high column density of 9 x 10^° cm^^ to match the minor axis emission profile. Comparison 
of modeled and data pv diagrams for other values of Rq constrain this parameter to be in 
the range 160-200" (22.8-28.5 kpc). Even for such a narrow ring, the high-latitude velocity 
profiles tend to be broader than in the data. For values of dVrot/dz outside the range -3 
to -4 km s~^ arcsec^^ the gradient in the halo component in the pv diagrams is not well 
matched. However, a somewhat higher gradient is required for the SE quadrant, and models 
indicate a value of dVrot/dz — —5 km s~^ arcsec"^. All halo models employ i — 86.25° to 
match the minor-axis emission profile. 

The narrowness of the ring also introduces a strong edge-brightening in the halo in 
zeroth-moment maps (Figure 12) due to the projection of the ring which is clearly not seen 
in the data. This halo would therefore need to be much more asymmetric than the warp 
of Models A and B. The normahzation of the modeled halo also now depends on the radial 
range used to form the minor axis emission profiles. We have chosen to normalize using the 
same radial range as above, realizing that the edge-brightening will lead to excess mass in 
the modeled halo. 
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An annulus wider than 1' (8.6 kpc), although featuring a more plausible column density, 
causes the high-latitude velocity profiles to broaden, if only slightly. Model E, where the 
annulus is 1.25' wide and its column density is 1.8 x 10^° cm~^, shows this effect. 

The major-axis pv diagram for the data and Models B and D are compared in Figure 
7. The general shape of the flat part of the pv diagram is reproduced quite well, given 
that we do not attempt to model the aforementioned asymmetry in the emission around the 
major axis. The general slope of the rising part of the diagram is underestimated. In the 
models, the slope is dictated by the fact that most of the gas is in the outer (i? = 3') ring 
in our axisymmctric radial profile. Asymmetries, in the sense of gas being more centrally 
concentrated in the vicinity of the inner ring at major axis offsets within 1' of the center, 
could be the cause of this discrepancy. 

The velocity dispersion is not well constrained because of the somewhat low velocity 
resolution. A value of 10 km s~^ was used in the models for both disk and halo. 

We roughly estimate the fraction of mass in the extended component in two different 
ways. First, we simply assume that in the zeroth-moment map of Figure 2, all emission 
above a certain height can be attributed to the extended component. Using two values of 
this height - 9" and 11" (1.3 and 1.6 kpc) -, its mass is 1.2 — 1.6 x 10^ Mq. The second 
estimate comes from the models. In the warp model, the excess column density in the 
warped region accounts for 1.8 x 10^ Mq. For the halo models, we calculate the mass in the 
halo component and find 3.1 x 10^ Mq. These are both overestimates, especially the latter, 
because the extra column density along the line of sight required by the axisymmetric models 
overestimates the emission along the major axis. The estimate from the data themselves is 
likely to be superior. 

Although both models are inaccurate because of their symmetry, the problem is much 
worse for the halo model. This model is also much more contrived in that it requires a very 
large mass in a very small radial range. The narrowness of the high-latitude spectra is also 
better reproduced by the warp model. We therefore conclude that the extended component 
is a warp rather than a halo. Further fine tuning of the warp model may produce a slightly 
better match to the data, but the asymmetries preclude a significant improvement, and hence 
we do not pursue this further here. 

3.2. Radial Inflow 

Figure 13 (top panel) shows a pv diagram along the minor axis of the galaxy. A shift 
in mean velocity between the positive and negative sides is apparent in the disk component. 



while a very small shift also exists in the warp component (the double-peaked nature of the 
disk component is due to the inner ring of 1.5' radius as discussed above). Such a shift is 
a potential sign of radial inflow in a not quite edge-on galaxy. All diagrams for major axis 
offsets within one beam width of this one show this signature. For the warp, such a shift is 
also seen, albeit at a slightly smaller amplitude, in the corresponding pv diagram of Model B 
(bottom panel of Figure 13), and we cannot rule out that it is an effect of the warp geometry. 
For the disk, since there is strong kinematic evidence of a bar (see §1), the signature would 
suggest radial inflow in the inner ring-like gas distribution surrounding the bar. To measure 
the inflow amplitude, we have created two spectra from the minor axis pv diagram, averaged 
over the range 3-18" in the minor axis direction on the east and west sides. A Gaussian is 
fit to each spectrum to find the mean velocity, and half of the velocity difference between 
the east and west sides is taken as the infiow velocity. The result is Vinfiow — 9.7 ± 0.4 km 
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Fig. 13. — Minor-axis position-velocity diagram for the full-resolution cube (top) and Model 
B (bottom). Contour levels are 1.5, 3, 6 and 12 times the 1 a noise in column density units 
of 4.6 X cm-2. 
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3.3. A Possible Supershell 

A structure in tlie SW quadrant above the disk has the form of a closed shell - better 
seen in Figure 14 - although it is always possible that projection of unrelated features has 
given it this appearance. The total mass is about 10® Mq, and the extent parallel to the 
major axis is about 3.4 kpc. Position-velocity diagrams parallel to the minor axis for the 
northern wall, center and southern wall of the shell are shown in Figure 15. Above about 
10" (1.4 kpc) from the plane, the walls show emission at the velocities of the warp as well 
as a broad component extending to about the terminal velocity in the northern wall, and 
some evidence for line splitting in the southern wall. The northern wall is detected over 
nine velocity channels, or about 190 km s~^. The line splitting in the southern wall is five 
channels or about 100 km s~^. Given that the velocities in the warp component are due to 
the warp geometry, it would be incorrect to conclude that this line broadening/splitting is 
due to expansion. Rather, it may just indicate that the feature is extended along t he line 



of sigh t. There is no indication of a corresponding ionized feature in the Ha image of iRand 



fll996h 
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Fig. 14. — Close-up of the zeroth-moment map of NGC 5746 showing the shell-hke structure 
described in the text. 
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Fig. 15. — Position- velocity diagrams parallel to the minor axis from the full-resolution cube 
for the northern wall, center and southern wall of the possible supershell. Contour levels are 
2, 4, 8, 16 and 24 times the 1 a noise in column density units of 4.6 x 10^^ cm~^as in Figure 
1. 
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3.4. A Clumpy Warp? 

Figure 2 shows that the high- latitude emission is partially resolved into individual 
clumps. Almost thirty such features have been identified by visual inspection and are listed 
in Table 2. This selection is not meant to be complete or statistically well defined, but is 
simply meant to give a first-order idea of the masses of well detected features. These range 
from about 4 x 10^ M© to 10^ M©, and their masses sum to 7 x 10^ Mq, or about half of the 
total mass of the high-latitude emission. The most massive, numbered 18, is the shell-like 
structure discussed above. Most are unresolved but a few have a filamentary appearance 
and lengths of 3-6 kpc. The line widths are difficult to estimate for many features because 
of the velocity resolution of 42 km s~^ (FWHM) and the low signal-to-noise ratio of some of 
the detections. For example, the significant emission for three of the best detected clouds, 
numbered 5, 6 and 8 in Table 2, is confined to two channels and therefore the lines are 
unresolved. Others, such as 11 and 28, show emission over 6 — 7 channels, or 125 — 145 km 
s~^, and it is not clear if these represent actual internal motions or reflect the extent of the 
features along the fine of sight, for example. 

The 60" -resolution zeroth-moment map (Figure 4) reveals two faint vertical extensions 
from the plane in the SE quadrant, reaching heights of 15 and 20 kpc. Summing their 
emission above where they appear to merge with the disk, we obtain masses of 7.9 x 10^ Mq 
and 6.8 x 10^ Mq, for the northern and southern feature, respectively. A faint extension is 
also seen off the NW end of the disk. Its mass is about 6.2 x 10^ Mq. The first-moment map 
shows that the SE features have about the same velocity as the part of the warp that they 
connect to with little change in velocity along their vertical extent. 

3.5. High Latitude Features Not in the Warp 

We have carefully combed the channel maps and pv diagrams to look for high-latitude 
emission which does not have the velocities of the warp component but which may not be 
evident in Figure 2. We have found four such features. Position- velocity diagrams for the 
first three, at the center of their extent along the major axis, are shown in Figure 16. The 
features' locations are also marked in Figure 3, and their properties are summarized in Table 
3. One is an unresolved feature in the NE quadrant of the halo. The second and third are 
diffuse features with appreciable extents along the major axis of 0.7' and 0.6' (6.0 and 5.1 
kpc), respectively. The second feature is also counterrotating and its emission spans more 
than 300 km s~^. The fourth is the low-T^e/ side of the possible shell (Cloud 18 in Table 2). 
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Fig. 16. — Position- velocity diagrams parallel to the minor axis from the full-resolution cube 
showing emission (indicated by arrows) from three high-latitude features that do not have 
the velocities of the warp. From top to bottom, these are Clouds 1, 2, and 3 in Table 3. 
Contour levels are 1.5, 3, 6, 12, 24, 36 and 48 times the 1 a noise in column density units of 
4.6 X 10^9 cm-2. 
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3.6. NGC 5740 



As this barred (jErwinll2005l ) Sb galaxy is beyond the half-power radius of the primary 
beam, the detailed structure may not be accurate and we limit ourselves to general properties 
only. We assume its distance to be the same as that of NGC 5746. Channel maps of the 
full- resolution cube are shown in Figure 17. Zeroth- and first-moment maps of NGC 5740 
are presented in Figure 18, using all data above 2a in two consecutive channels. The zeroth- 
moment map is overlaid on a Digitized Sky Survey red image in Figure 19. There is a disk 
of about 30 kpc diameter and a large extension to the NNW. The disk has a mass of about 
5.3 X 10^ M0 and the extension 9.4 x 10® M©. For comparison. ISchneider et al.l (119861 ) quote 
a mass in atomic hydrogen of 7.0 x 10^ Mq, scaled to our distance. 
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Fig. 17. — Channel maps from the full-resolution cube for NGC 5740, corrected for primary 
beam attenuation. Contour levels are 3, 6, 12, 24, 36 and 48 times 4.6 x 10^^ cm~^. Each 
map is labeled with its heliocentric velocity. 
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Fig. 18. — Left: Zeroth-moment map of NGC 5740 from the full-resolution, primary-beam- 
corrected cube. Contour levels are 2, 4, 8, 12 and 16 times 10^° cm~^. Right: First-moment 
map of NGC 5740. Contours run from 1440 to 1740 (the closed contour on the NE side at 
Decl. 1° 42') km s"^ in 20 km s~^ increments from SW to NE. 
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Fig. 19. — Zeroth-moment map of NGC 5740 from the full-resolution cube overlaid on a 
Digitized Sky Survey red image. Contour levels are as in Figure 18 (left). 
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A ROTCUR analysis of the disk indicates a PA of 339 ± 0.5°, an inclination of 58 ± 2°, 
Vsys = 1586 ± 6 km s-\ a dynamical center of R.A. 14^^ 44°^ 24.4^ Decl. 1° 40' 52", with 
approximate uncertainty 3" in R.A. and 8" in Decl. The first-moment map suggests regular 
rotation continues into the NNW extension. Keeping Vsys and the dynamical center fixed, 
and fitting the receding half of the galaxy only, the extension exhibits little change in PA 
but a continued fall in Vrot and an increase in inclination to 61° at 2' and 70° at 4'. The 
derived rotation curve is shown in Figure 20. The dynamical mass within R = 2' (17.1 kpc) 
is 1.2 X 10^^ M0. 

Morphologically, the NNW extension is similar to extensions in s everal Virgo Cluste r 



galaxies, most of which are best explained by ram pressure stripping (IChung et al.l 120071 ). 



There is no obvious evidence of compression of the gas on the o pposite side o f the d isk, 



however. The condition for ram pressure stripping, as given by iVoUmer et al.l (120011 ) . is 
PiGMVgg^i > ^isAiVrot/^y where pigm is the intragroup medium (IGM) density, and Vgai is 
the velocity of the galaxy through the medium, for gas at surface density T,jsm at a radius 
R from the galactic center. At the outskirts of the main disk the typical surface density 
is about 2 MqPc"^, while Vrot = 175 km s~^ The condition for stripping is then (n/10~'^ 
cm~^)(f/100 km s~^)^ > 14. Both required numbers are not constrained. We find that the 
one-dimensional velocity dispersion for the catalogued galaxies of the NGC 5746 group is 120 
km s~^. The IGM density is unknown . X-ray bright groups are found to be Hl-deficient by 



Sengupta fc BalasubramanyamI (120061 ). who calculate a range of IGM densities from 5 x 10~ 
to 2 X 10^^ cm^^, but it is unknown how extensive X-ray emission is in the NGC 5746 group. 
Hence, it is unclear whether the stripping condition can be met. Tidal stripping seems less 
likely given the morphology and that the closest group galaxy on the sky, NGC 5746, is 155 
kpc away in projection. 
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Fig. 20. — Filled circles show the rotation curve for NGC 5740 derived from a tilted ring 
analysis. Open circles show the extension of the rotation curve into the NNW extended 
emission component. 
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Discussion 



Almost all of the high-latitude gas we have detected can be associated with the warp. 
At the level of sensitivity of our data, only the features listed in Table 3, summing to about 
10® Mq, may have originated in a disk-halo flow or infall. The heights off the plane of their 
centers range from 2 to 8 kpc. One of the more massive ones is counter-rotating and is 
unlikely to have originated in the disk. A few x 10^ Mp^ worth of counter-ro tating clouds 
were also discovered in the neutral halo of NGC 891 by lOosterloo et al.l (120071 ). 



Are these few clouds consistent with predict ions from the ha l o ther mal instability sim- 
ulations mentioned in §1? In the simulations of ISommer-LarsenI ( l2006l ). most of the warm 
clouds are within 50 kpc from the center of the galaxy and total 10® M0, with most having 
masses a few times 10^ to a few times 10^ Mr:^ for a galaxy like the Milky Way. By compari- 
son, for an M33-like galaxy, iKaufmann et al.l (120061 ) expect clouds of radii 0.1 — 0.6 kpc and 
masses from 10^ to a few times 10^ M©, totaling ~ 10^ — 10® Mq depending on the simulation 
parameters, confined to within about 10-20 kpc above the disk. If the few clouds in Table 3 
are due to infall, we note that their total mass is comparable to these model predictions, with 
the caveat that the fraction of neutral gas is uncertain in the simulations. The small number 
of detected clouds is i n fact comparable to that in the "very high resolution" simulation of 
Sommer-Larseru (120061 ). where the mass resolution is comparable to that of our data. 



For galaxies like the Milky Way, Mailer fc BuUockl (|200J) predict 2x10^° Mq of mostly 
ionized clouds of mass 5 x 10^ Mq and size ~ 1 kpc, extending to 150 kpc from the center 
of the galaxy. At our sensitivity, it is difficult to rule out with confidence such a widespread 
population of clouds from the HI data alone. 

It would be interesting to know how much mass in warm clouds forms in simulations 
of more massive systems like NGC 5746, and where it is located. The cooling rate, at least, 
of the hot gas in this galaxy can be est i mated from the observations, but is very uncertain. 
The rate inferred by iRasmussen et al.l (120061 ). who exclude the projected disk area where 
there may be contamination by other sources of hot gas, is only 0.2 Mq (yr)"^. The rate is 
modest because the estimated halo temperature is rather high, at 0.56 keV. However, most 
of the halo cooling is likely to be occurring in the excluded area, and their models of massive 
disks suggest that the cooling rate could be 5-10 times higher as a result. In summary, it 
is unclear whether the observed clouds are consistent with theoretical expectations for halo 
cooling in a massive disk galaxy like NGC 5746. 

What fraction of these halo clouds do we expect to be fully or mostly ionized? This 
depends on the incident ionizing fiux, the mass, and most importantly, the bounding pressure 
on the cloud. If we assume that the ionized gas in a cloud is in thermal pressure equilibrium 
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with an outside pressure, the density of a cloud is given by n_2 = 4:.55T^^ [ps/k], where 
n_2 = n/^/10~^ cm~^ is the normahzed hydrogen density of the cloud, T4 = T/(10^ K) 
is the normalized temperature, and [ps/fc] = [p//c]/(10^cm~^ K) is the normalized thermal 
pressure. (We use here p/rin = 2.4 x 10~^^ for solar metallicity gas, ritot = 2.2?t,j:/, and 
Ue = 1.2nH for fully ionized plasma). For a constant density spherical cloud, the cloud 



radius is given by i? = (532 pc)Mq^^tI^^ [ps/k] 



-1/3 



We assume that the outer skin of this cloud will be ionized by a (plane-parallel) flux 



(f)/ (10^ pho tons cm 



The metagalactic flux is estimated 



of H ionizing pho tons, (t>^ 
to be 05 < 0.45 (ISternberg et al.l |2002| . and references th erein), while high velocity clouds 
in the Galactic halo experience 05 = 2 JXufte et al.ll20oi ). The thickness of this skin, £, is 



determined by the balance of ionizations and recombinations, = a^rigf where the Case B 
recombination coefficient is = 3 x 10^^^ cm'^ s~^. (Note that for optically thick clouds, 
(f) = as EM, so that the surface averaged emission measure is a simple measure of the 
ionizing flux.) In this case, the skin depth is given hj i = (35.7 pc)05T4 [ps/k]~'^. 

With these approximations, the ionized volume and mass of the cloud is given by fion = 
1 - [1 - {i/R)f or fion = 3{£/R) - 3{i/RY + ... for small (i/R). The condition for a cloud 
to be fully ionized is 



< (197 cm~^ K) (Pl^^TiM'''-". (1) 
For very low ionization fraction {fion < 5%), this criterion is well approximated by 



'1/3 



(380 cm-3 K)/, 



-3/5,3/5 



TaM, 



-1/3 



so that an ionization fraction of less than 1% requires [p/k] > (6022 cm ^ K) 



(2) 



The se nsitivity of the io i iizatio n fraction of ga s cloud s to t he outside pressure h a s been 
studied by iFerrara fc Fkldl Jl994h . IWolfire et aP Jl995h . and iMalonev fc PutmanI J2003h 
with comparable results. The normalization constant relating the pressure to the ionization 
fraction and cloud mass will depend on the assumed geometry. The presence of central 
condensation due to dark matter halos ca n also lower somewhat the pressure necessary to 
form neutral clouds (ISternberg et al.ll2002l ). 



These considerations show that not only are X-ray halos associated with galaxy forma- 
tion more likely to form condensations, but that the higher pressures of such halos result in 
a much larger neutral fraction for the c ondensation. Using the density and temperature esti- 
mate for the X-ray halo of NGC 5746 fiRasmussen et al.ll2006l ). the halo pressure is given by 
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Ip/k] = 6200 cm"^ K. Such a high bounding pressure would guarantee that clouds embedded 
in this hot gas will be primarily neutral. It is also interesting to note that the bounding 
pressure for the cloud can be either the thermal gaseous halo pressure or the ram pressure 
associated with a cloud's motion through the halo. The one cloud showing counter rotation 
presumably has a much higher ram pressure than most of the rest of neutral gas detected 
and could be expected to have a higher neutral fraction. 

At larger radius than that traced by the X-ray gas, it remains possible that the pressure 
will drop to such a point that a large population of fully ionized clouds might be present, 
but uncertainties in the X-ray data and analysis prevent such pressure variations from being 
assessed. Although the predicted emission measure from the front and back face of an 
individual cloud is a relatively low EM= 2 x (0.11 pc cm~®) 05, a population of unresolved 
smaller clouds in the beam might boost the Ha emission to detectable levels. However, 
processes such as thermal conduction and eyapora tion should set a lower limit on the mass 



of clo u d tha t can survive (IMaller fc Bulloclal2004l ). We have reexamined the Ha image of 



Randl (119961 ). smoothed to 1.5" resolution - yielding a noise level of 2.2 pc cm ^ - and find 



no evidence for such bright clouds in the halo. 

To attempt to constrain the origin of the neutral gas halo in this and other galaxies, 
we can compare the neutral halo mass and mass fraction with estimates for other galaxies 
and look for correlations with other parameters (Table 4). In this comparison one should 
keep in mind that some of the galaxies are edge-on, one (NGC 6946) very much face-on, 
and some at intermediate inclinations (NGC 4559, NGC 2403 and NGC 253), and thus the 
characterization of halo gas necessarily varies from galaxy to galaxy. In the galaxies not 
viewed edge-on, the extraplanar gas is identified by its anomalous velocities compared to 
the bulk of the emission. We also note that the sensitivity and linear resolution of these 
observations vary significantly. 

NGC 5746 ranks near the bottom of the table in terms of halo mass and mass fraction 
and also has one of the lowest levels of star formation activity, as estimated from Lpm/ D^^, 
lending support for an origin in a weak disk-halo flow. However, there is no discernable trend 
among the galaxies in this direction. NGC 6946 and NGC 253 stand out as having little HI in 
their halos given their star forming activity. The halo HI in the latter is part icularly puzzling 



as it is found only on one side of the disk, away from the center of the galaxy (IBoomsma et al. 



20051 ). It is possible that the neutral gas has been swept radially outwa rd by the c entra l 



starburst outflow or that much of the halo gas is ionized. In NGC 6946, iBoomsmal (120071 ) 
note that the mass of halo gas more than doubles in a map at 64" resolution vs. one at 22" 
resolution. The rotation speeds of the galaxies span a factor of three, and thus the ability of 
supernova power to raise gas off the plane should vary significantly from galaxy to galaxy. 
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In particular, the high fractions of halo HI in NGC 4559, NGC 2403 and UGC 7321 could 
be in part due to their low mass. Finally, there does not seem to be any correlation between 
halo mass or mass fraction and rotation speed alone. 

However, at this point it is not unreasonable to conclude for NGC 5746 that the three 
non-counterrotating clouds in Table 3 may just as likely be due to a relatively inactive disk- 
halo cycle as to infall. One possible further clue to their origin would be infrared emission, 
and indeed imaging observations of this and other edge-ons with the Spitzer Space Telescope 
have been proposed. A detection of dust would argue against an origin in low-metallicity 
inf ailing gas. 

A larger sample of edge-ons needs to be observed with comparable sensitivity and angu- 
lar resolution to understand how the mass and kinematics (rotational lags in particular) of 
neutral halo gas relates to other galactic properties such as star formation activity, assuming 
that such gas can always be kinematically distinguished from warps and flares. A better 
measure of such activity than the somewhat crude tracer e mployed here would be the 24/im 
surface density as measured by the Spitzer Space Telescope (ICalzetti et al.ll2005l ). hence 24/xm 
maps of a large sample of edge-ons would be of great benefit in such comparisons. For NGC 
5746, deeper observations with higher spectral resolution are clearly called for in order to 
carry out a more complete census of halo gas that can be distinguished from the warp, so 
that its spatial distribution, kinematics, and cloud masses can be compared to the various 
models discussed here. 



It should also be remembered that NGC 5746 is in a group environment, and some of 
the high-latitude gas may owe its origin to previous encounters, although the apparently low 
metallicity of the hot gas argues against encounters with other large galaxies. 

We reemphasize that our best warp model (Model B) is not necessarily a unique fit to 
the data, especially because of asymmetries, but does demonstrate that a warp is a superior 
explanation to a lagging halo. It is not surprising that a wa rp is present in NGC 5746, as they 



are common in disk galaxies (e.g iGarcia-Ruiz et al.ll2002l ). As far a s the warp p arameters 



are concerned, we note that it is generally the case in such galaxies (lBriggslll990l ) that the 
line of nodes increases in the direct ion of galaxy ro tation for radii beyond the Holmberg 
radius, which is 4.5' for NGC 5746 (jHolmberglll958r). Our warp mode l implies an increase 



beginning just inside this radius. More recently, Ivan der Kruit 



(120071 ) finds evidence from 



edge-ons that warps begin just beyond the radius where a truncation is evident in the stellar 
disk, but finds no indication of such a truncation in NGC 5746. 



Finally, there is the possibility that the warp itself is caused by infall. The origin 
of warps is debated, but the possibility receiving most attention recently is that they 
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are due t o infalling gas w i th an angular momentum ve ctor tilted with respect to the in- 
ner disk (Jiang fc BinnevI Il999l : IShen &: SellwoodI |2006| ). In the N-body simulations of 



Shen &: SellwoodI (120061 ) . in which particles are injected into an outer torus inclined at a 



fixed angle with respect to the main disk, the warp can persist for several Gyr, making it 
difficult to infer the accretion history from the mere presence of a warp. That is, the large 
gas mass in the warp may refiect accretion over a substantial period of time and not nec- 
essarily demand a high accretion rate, as would more likely be the case such a large mass 
were contained in a lagging halo. On the other hand, if the centrally concentrated X-ray 
halo emission reflects the hot gas density, then one would expect any warm clouds recently 
formed by cooling out of the hot phase to be more centrally concentrated. 



5. Conclusions 

We have observed the massive edge-on spiral NGC 5746 with the VLA, in an attempt to 
find a vertically extended component of atomic gas as predicted by recent galaxy formation 
models. A high-latitude component has been discovered, but almost all of its mass, summing 
to 1.2 — 1.6 X 10^ Mq - or about 15% of the total HI mass - is more readily explained as 
a warp than as a halo. The warp must be asymmetric in its column density distribution 
and geometry, but this is not unusual. If the high-latitude component is a lagging halo, it 
must feature a large mass in a rather narrow radial range, be centered at a large radius, 
and have a much larger degree of asymmetry than the warp. Even then, such a model has 
difficulty reproducing the narrowness of the high-latitude spectra as well as the warp model. 
We therefore conclude that a warp is the more likely explanation. The warp itself may be a 
result of infall according to recent models. It partially resolves into clouds of mass 4 x 10^ 
Mq to 10^ Mq, accounting for about half the mass in the warp. 

We have found four high-latitude features at velocities distinct from the warp, totaling 
about 10® Mq. These could be accreting onto the disk, or they may originate in a disk- 
halo fiow, although there is no other sign of such a fiow in this galaxy with little star 
forming activity. One cloud is counter-rotating and must have an external origin. The 
cl oud properties are rough ly c omparable to those exp ected in the galaxy formation models 



of iKaufmann et al.l (120061 ) and ISommer-LarsenI (120061 ) where warm clouds form by thermal 



instabilities in the residual hot halos. However, these models have not been tailored to 
galaxies as massive as NGC 5746, and it would be interesting to compare our results with 
such a model. Given the calculated pressure in the hot halo, we expect any clouds to be 
primarily neutral. Alternatively, we cannot rule out a weak disk-halo fiow as the origin for 
these clouds, except for the counter-rotating one. How neutral gaseous halos relate to other 
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galaxian properties needs to be tested with more observations. 

The disk shows the signature of radial inflow at the level of about 10 km s~^. This 
is most likely due to the bar, which is known to have a strong kinematical signature from 
optical long-sht spectra. 

The group member NGC 5740 is also detected, albeit outside the half-power points of 
the primary beam. The most important result here is a broad extension of gas to the NE, 
which is more likely to be due to ram pressure than tidal stripping, although whether the 
condition for such stripping is met is unknown. 

We thank J. Mulchaey for discussions regarding X-rays in groups, and L. Sparke for 
discussion about warps. We also thank an anonymous referee for many useful comments. 
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Table 1. Models of the High Latitude Emission 



Model 


Description 


Incl. 


P.A. 


Disk rotation 


High-latitude 






(degrees) 


(degrees) 


curve 


Nh profile 


A 


warp along 
line of sight 


Figure 11 


350 


Figure 8 


Figure 6 


B 


warp at angle to 
line of sight 


Figure 11 


Figure 11 


Figure 8 


Figure 6 


C 


lagging halo 
(radial profile of disk) 


86.25 


350 


Figure 8 


Figure 6 


D 


lagging halo 
(narrow ring) 


86.25 


350 


Figure 8 


Figure 6 


E 


lagging halo 


86.25 


350 


Figure 8 


Figure 6 



(broader ring) 
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Table 2. Properties of Clumps in Extended Component 
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'^Numbers generally run in a counter-clockwise direction around the disk starting at the northern end. 
''Offsets are from the dynamical center determined in the text. 





Table 3. 


High- 2; 


Emission Features Not Associated with the Warp 




Number 


R.A. 


Dccl. 


Offset along major axis" 


Offset along minor axis'' 


Mass 


Range in Vhel 




(J2000.0) 


(J2000.0) 


(arcmin) 


(arcmin) 
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(km s^-"-) 
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-0.70 
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1391-1495 



* Offsets are from the dynamical center determined in the text. 
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Table 4. Neutral Gaseous Halos of Nearby Galaxies 



Galaxy HI halo mass Percentage of HI Lpir/D2^ Rotation speed^ 

(10^ Mq) mass in halo (10'"^ erg s"! kpc'^) (km s"!) 



NGC 891'' 


12 


30 


3.5 


225 


NGC 4559^^ 


6 


10 


0.8 


120 


NGC 3044'* 


4 


7 


4.7 


150 


NGC 2403== 


3 


10 


1.1 


135 


NGC 6946* 


3 


4 


9.6 


170 


UGC 7321S 


1.3 


12 


0.1 


105 


NGC 5746 


1.0 


1 


0.4 


310 


NGC 253'^ 


0.8 


3 


10.6 


230 


M3l' 


0.35 


1 


0.2 


270 



^An approximate value for the flat part of the disk rotation curve. 
''Halo mass and rotation speed from lOosterloo et al.l (|2007| ) 
'^Halo mass and rotation speed from lBarbieri et al. I ll2005h 
''Halo mass and rotation speed from lLee fc IrwinI l ll997l) 
"Halo mass and rotation speed from iFraternali et al.l (|2002| ) 
'Halo mass and rotation speed from iBoomsmal ||2007| ) 

^Halo mass from lMatthews &: WoodI l l2003h and rotation speed from lUson fc Matthew3 | |2003|) 
''Halo mass and rotation speed from lBoomsma et al. 1 ll2005h 

'Halo mass from iThilker et al.l |2004l '). total mass from iBraun fc Walterbo i l ll992h . rotation 
speed from lBrinks fc BurtonI lll984 ) 



